Azimuthal dependence of collective expansion for symmetric heavy ion 

collisions 
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Detailed studies of the azimuthal dependence of the 
mean fragment and flow energies in the Au+Au and 
Xe+Csl systems are reported as a function of incident 
energy and centrality. Comparisons between data and 
model calculations show that the flow energy values 
along different azimuthal directions could be viewed as 
snapshots of the flreball expansion with different expo- 
sure times. For the same number of participating nu- 
cleons more transversally elongated participant shapes 
from the heavier system produce less collective trans- 
verse energy. Good agreement with BUU calculations is 
obtained for a soft nuclear equation of state. 
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One of the main motivations to study heavy ion 
collisions at high energy is to obtain information on 
the equation of state (EoS) for nuclear matter un- 
der conditions of pressure and temperature differ- 
ent from those in normal nuclei. The search for hot 
and dense nuclear matter created in such collisions is 
confronted with dynamical consequences of the high 
incident energy necessary to reach such conditions 
and with the difficulty to reach a thermal equilib- 
rium in finite systems. Dynamical aspects refer not 
only to the initial phase of the collision, but also 
to the evolution stage of the formed fireball. Thus, 
detailed experimental information on the expansion 
dynamics is required. The simplest situation cor- 
responds to central collisions with the advantage of 
the azimuthal symmetry and of the lack of specta- 
tor matter. Predicted in early seventies [1,2], the 
collective expansion of the hot and dense fireball 
produced in central collisions was evidenced experi- 



mentally [3-11]. Although central collisions seem to 
deliver the cleanest signal on the collective expan- 
sion on first sight, two issues are worth mentioning: 
i) While the axial symmetry of the dynamical evo- 
lution holds, the spherical symmetry has to be in- 
spected. Preequilibrium emission and transparency 
effects could influence the spherical symmetry of the 
expanding system, ii) With regard to reaching pres- 
sure, the nuclear matter, not being confined in trans- 
verse directions, can escape freely in any direction 
perpendicular to the collision axis starting from the 
very first moments of the collision. For reduced cen- 
trality, other complications appear. One has to deal 
with rotating expanding objects in the presence of 
spectator matter. Nevertheless there are also some 
advantages in studying less central collisions: i) ro- 
tation and shadowing can be used as internal clocks 
for getting deeper information on the expansion dy- 
namics, ii) the centrality can be used to control the 
shape and content of the fireball and of the shadow- 
ing matter, iii) for a given centrality the passage time 
of the shadowing objects can be controlled varying 
the incident energy, iv) the confinement of the spec- 
tators becomes more efficient in transverse direc- 
tions than in the central collisions. Symmetry con- 
siderations imply two dominant components in the 
particle transverse emission: azimuthally symmetric 
emission and an elliptic modulation of that emission 
(squeeze-out) which has been predicted by hydro- 
dynamical calculations [12] and extensively studied 
experimentally [13-24]. The squeeze-out has been 
studied, in particular, as a function of centrality, 
type of emitted particle, transverse momentum and 



mass of the colliding systems. A considerable collec- 
tivity was identified in the transverse emission pat- 
tern and it was found, from fitting the energy spectra 
with a radially expanding source, that the collectiv- 
ity itself exhibits an elliptic modulation [19]. The 
present paper is devoted to the elliptic modulation 
of collectivity as reflected in the dependence of av- 
erage fragment energies on the fragment mass [22] . 

Notably, without any modulation of collectivity, 
the squeeze-out itself could just represent pure geo- 
metric shadowing of particle emission from the par- 
ticipant zone. Significant modulations, as a counter- 
part, should show significant influence by early pres- 
sure, with the compression and excitation energy 
getting converted into collective energy at the stage 
of high density when the spectators are present. 
Comparisons of different centralities and system 
masses arc needed in order to imderstand how the 
geometry of the participant and spectator zones af- 
fects the collective expansion. Excitation functions 
for the collectivity anisotropy are mandatory to un- 
derstand how the collective expansion builds up at 
different densities and excitations. 

The present work is based on ^^^Au +^^^Au and 
i^OXe -I- I33csi27i data from FOPI Phase II exper- 
iments at the SIS of GSI. Detector details can be 
found elsewhere [25,26,24]. The main component 
implemented in the FOPI configuration is a Central 
Drift Chamber (CDC) [26]. Two criteria for central- 
ity selection were combined: the particle multiplic- 
ity, with a higher selectivity at large impact parame- 
ters and the ratio of transverse and longitudinal en- 
ergies in the cm. system Erat = E±,t/ J2i ^^i, 
as a better filter for more central collisions. Two 
windows in the CDC multipUcity, CM2 and CMS, 
have been used to select impact parameters in the 
range of 6-8 fm and 4-6 fm. Similarly, the ER4 and 
ER5 windows in Erat select the ranges 2-4 fm and 0- 
2 fm, respectively. The impact parameter estimates 
and the estimates of the number of participating nu- 
cleons Apart are based on the geometrical "sharp 
cut-00" approximation. For reaction plane determi- 
nation, the transverse-momentum analysis method 
has been used [27]. Studies of the squeeze-out phe- 
nomena revealed from the very beginning [13] the 
importance of performing the azimuthal distribu- 
tion analysis in a reference frame with the polar axis 
along the sidewards flow direction. The present anal- 
ysis used a reference frame in which the ellipsoidal 
pattern of the azimuthal particle distribution at mid- 
rapidity maximizes the ratio of the two transverse 
semiaxes [19,22,24]. The previous comprehensive de- 
scription of the elliptic modulation of collectivity in 
the beam range 0.25 - 1.15 A-GeV [19] was based on 
a parameterization of the deuteron (A=2) tranverse- 
mass spectra with an expression characteristic of ra- 
dially symmetric shell expansion. The use of such 
an expression for studying the modulations is inher- 
ently contradictory and implies a specific model. We 



prefer to present experimental information free from 
any model. 

For the mentioned reasons we concentrate on the 
mean kinetic energy in the cm. system, (E^™), and 
on the flow energy extracted from its dependence 
on the mass of the reaction products (Z=l,2 and 3) 
within a polar angular range of 80° < 9cm <100° [22]. 
In order to extract {E^i^), one needs complete en- 
ergy spectra. Due to the peculiar shape of the shad- 
ows of subdetector borders in the rotated reference 
frame, the experimental spectra have been analyzed 
in the azimuthal ranges [0°, 90°] and [270°, 360°]. 
These two ranges have been overlapped to decrease 
statistical errors. They are plotted as full symbols 
in Fig. 1 for five bins from 0° to 90° and were then 
reflected (open symbols) in order to generate the full 
angular range of 0°-360°. 
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FIG. 1. Mean cm. energy (S^^) of Z=2 products as 
a function of the azimuthal angle, shown for: Au-|-Au 
at ER4 ccntrality and different energies (left); Au-fAu 
at 250 A-McV and different centralities (middle), and 
Au+Au and Xe+Csl at 250 A-McV and ER4 centrality 
(right). The error bars include systematic effects. 

Figure 1 presents, as an example, the azimuthal 
dependence of {E^^) for Z=2 products as a function 
of the incident energy in Au -|- Au at ER4 central- 
ity, as a function of centrality in Au -|- Au at 250 
A-MeV, and for the two measured systems at 250 
A-MeV and ER4 centrality. The presented informa- 
tion is independent of the anisotropy of the yield dis- 
tribution. The main contribution to the error bars 
comes from systematic effects, the statistical ones 
being at the level of symbol sizes. We had to com- 
bine information from two subdetector systems and 
this was done by looking at all fragments as a func- 
tion of their charge. In doing this, we had to correct 
the Z=l and Z=2 energy spectra from the CDC, 
based on previous FOPI data [8], in order to take 
into account the fact that the ^He fragments are 
not well separated from the tritons in some parts of 
the momentum space. The azimuthal asymmetry of 
i^kiii) increases as a function of incident energy and 
mass of the colliding nuclei, reaching a maximum in 
mid-central collisions. The value averaged over the 
azimuth increases with the beam energy, centrality 
and mass of interacting system. As shown by the 



solid lines, the data follow very well the behavior 
(^fem) = Eli^ - ^Ek„r-cos2^. Fits to the depen- 
dence of mean kinetic energies (i?^™) on fragment 
mass A with the nonrelativistic expression 
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yield the average flow energy per nucleon Ecou and 
the "temperature" "T" (mo is the nucleon rest 
mass). Lack of an explicit treatment of the Coulomb 
effects leads to a systematically overestimated value 
of the real temperature T [8] . 
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FIG. 2. E°„u (top-left panel), E^Ji" (top-right panel) 
and AEcoii (bottom) as a function of Apart, correspond- 
ing to the four centralities mentioned in the text, at dif- 
ferent incident energies in the Au-|-Au (full symbols) and 
Xe+Csl (open symbols) systems. Second-order polyno- 
mial fits, represented by dashed lines, serve to guide the 
eye. 

The flow energy E,,oU=l/2-mo< > and "T" 

can be fit nicely with EcoU = E^^^ - AEcoU ■ cos 2$ 
and "T" = "To" - A"T"-cos2$, respectively. Ecoii 
exhibits a strong elliptic anisotropy, with the largest 
energy values in the direction perpendicular to the 
reaction plane. Both the average over the azimuth of 
the collective energy E^^^ and the elliptic anisotropy 
AEcoii increase continuously with the incident en- 
ergy over the studied region. The temperature pa- 
rameter " T" stays roughly constant as a function of 
the azimuth at three lowest beam energies and de- 
velops oscillations A"T" in the range 1-2.5 MeV, at 
250 A-MeV. These could be indicative of the small 
variations of temperatures at the different sides of 
the participant fireball [7]. Figure 2 provides global 
results: E%, E-^o?fiEcoU at 90°) and AE.oU as 
functions of Apart for both measured systems. At 
90 A-MeV, the in-plane and out-of-plane flow val- 
ues are very similar (AEcoii^O) at all centralities in 
Au-f Au, paralleling the observations for the yields 
for which a transition from in-plane to out-of-plane 
enhancement was found [24] as a function of inci- 
dent energy. At all energies, E^^^ increases with the 
centrality corresponding to an increasing baryonic 
number A„„,rt of the fireball. Although the error 



bars, which include the systematic effects, are large, 
the relative errors for the data at different Apart are 
small, of the order of 0.08 MeV. At higher energies, a 
maximum of AEcoii in mid-central Au-|-Au collisions 
becomes visible. A comparison of the two systems 
250 A-MeV in the 90° direction, outside of the spec- 
tator shadow. Fig. 2 top-right, shows a difference of 
the order of 20% in E^ff^ between the two systems, 
at the same Apart- One should notice that, for the 
same Apart, the fireball produced in the Xe -I- Csl 
case is more spherical than that in the Au -I- Au case. 
This remarkable dependence of the transverse fiow 
energy on the shape of the expanding nuclear zone 
is reported here for the first time. Flatter fireballs 
yield on the average less transverse expansion than 
more spherical ones. 




2 4 f, 8 10 12 14 1(> 

tJfm/c] 



12 _,;„14 16 

f" Jfm/c] 



§ 251 



Uj 



H0.0fm/c<t„<1 


6.5 fm/c 




A 14.0 fm/c<tt< 


6,5 fm/c 


■ 


T 15.5 fm/c<t.< 


6,5 fm/c 






■ 








T 




T 




! 








c) 




2 


4 


'A 

"IMF 



Uj 



4+ 



+\ /* ^ - 

d) 



50 100 150 200 250 300 350 



i^[deg] 

FIG. 3. a) Li-like fragment yield distribution in the 
plane of Ekin vs break-up time tb in the hybrid model, 
b) Mean kinetic energy of Li-like fragments emitted after 
a break-up time f™*". c) Mean kinetic energy of Z = 1, 2 
and 3 fragments as a function of their mass for different 
ranges of the break-up times, d) Collective energy as a 
function of azimuthal angle in the hybrid model. 

To get an insight on the main mechanism be- 
hind the observed experimental trends, we have 
calculated the expansion dynamics within a semi- 
analytical hybrid model [7] for a 200-nucleon fireball, 
which roughly corresponds to Apart at CMS central- 
ity in Au -I- Au collision. Notably, this model, in 
which expansion dynamics is combined with statis- 
tical features of cluster formation at freezeout, pre- 
dicts a decrease of the collective energy and of the 
temperature as the system deexcitation progresses. 
A two dimensional yield distribution as a function 
of Ekin and break-up time th for Li fragments in 250 
A-MeV Au+Au collision is presented in Fig. 3a. In 
the directions outside of the spectator shadow, i.e. 
perpendicular to the reaction plane, fragments with 
kinetic energies corresponding to all break-up times 
will be visible. The situation changes if the observer 
views the reaction from the reaction plane. Particles 



emitted directly from the fireball are visible once the 
spectators moved apart from the collision zone by a 
distance corresponding to at least half of the transi- 
tion time. Those particles which are emitted earlier 
will get redirected. At later times, the Ekin distri- 
butions become narrower following the decompres- 
sion and this effect gets more pronounced for heavier 
fragments. Figure 3b shows the mean energy of Li- 
like fragments as a function of a starting time 
for emission, from integrating over the yield in panel 
a) . Mean kinetic energies as a function of fragment 
mass number A, for different emission intervals in tb, 
can be seen in Fig. 3c. If one calculates the emission 
in a simple geometric picture, under the assumption 
that the centers of original nuclei pass each other in 
the middle of the expansion time shown in Fig. 3a, 
one gets the behavior of flow energies represented in 
Fig. 3d. Although quite simplistic, the model nicely 
reproduces the qualitative trends seen in the data, 
showing that the flow energy values could be viewed 
as snapshots of the fireball expansion dynamics with 
different exposure times for different azimuthal di- 
rections. 
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FIG. 4. -E^i„ and AEkin as a function of Apart, for 
Z=l (A=l,2,3) fragments, Au-fAu at 400 A-MeV. The 
experimental results are represented by triangles, while 
the BUU results are represented by gray zones for soft 
EoS and by dashed zones for stiff EoS, respectively. 

It is obvious that preequibrium processes can be 
important and that dynamic long-term anisotropics 
can be produced in the collective expansion with 
either effect being absent in the simple scenario 
above. Thus, a comparison with the ab ini- 
tio microscopic transport model becomes impor- 
tant. Correspondingly, for the 400 A-MeV Au-|-Au 
collisions, we present in Fig. 4 the results of 
the Boltzmann-Uehling-Uhlenbeck (BUU) trans- 
port code [28] using momentum dependent mean 
fields (m*/m=0.79), in-medium elastic cross sec- 
tions (ct — do tanh(cr-^'''^'^/(To) with ao — p~'^^^) 
and soft (K=210 MeV, gray zone) or stiff (K=380 
MeV, dashed zone) EoS. The light fragments (up to 
A=3) are produced in a few-nucleon processes in- 
verse to composite break-up. The measured relative 
yields are nicely reproduced, especially at higher in- 
cident energies [8]. The calculated yields have been 
smeared according to the measured reaction-plane 
dispersion. For (Ekin), little sensitivity to EoS is 
found, with either EoS parameterization yielding a 
quite good agreement with the data. As far as the 



values of AE^m are concerned, the calculations with 
the soft EoS reproduce the overall trend of the exper- 
iment while the calculations with the stiff EoS over- 
estimate significantly especially at higher centralities 
(lower impact parameter) the measured values. 

In summary, we presented results on the az- 
imuthal dependence of mean fragment and flow ener- 
gies in two symmetric systems, for different central- 
ities and incident energies. Corroborated by model 
estimates, indications emerged that different regions 
of the azimuth capture different periods of the cen- 
tral fireball expansion. In comparing results from 
two symmetric systems, it was possible to evidence 
that a more spherical fireball produces more trans- 
verse flow at a given participant number than a de- 
formed fireball. Comparisons with transport code 
predictions demonstrated that the soft EoS, in com- 
bination with a momentum-dependent mean fields 
and with in-medium cross sections, gives a good 
agreement with the experiment. 
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